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ABSTRACT 

We have obtained spectra of the J— 2-1 and J=10-9 transitions of cyanoacetylene (HC3N) toward 
a conection of positions in the most prominent filament, B213, in the Taurus molecular cloud. The 
analysis of the excitation conditions of these transitions reveals an average gas II2 volume density of 
(1.8±0.7) X 10^ cm~^. Based on column density derived from 2MASS and this volume density, the line 
of sight dimension of the high density portion of B213 is found to be ~ 0.12 pc, which is comparable 
to the smaller projected dimension and much smaller than the elongated dimension of B213 (~2.4 pc). 
B213 is thus likely a true cylinder-like filament rather than a sheet seen edge-on. The line width and 
velocity gradient seen in HC3N are also consistent with Taurus B213 being a self-gravitating filament 
in the early stage of either fragmentation and/or collapse. 

Subject headings: ISM:structure - ISM: individual objects:Taurus-ISM:molecules- 
Techniques : spectroscopic 



1. INTRODUCTION 

The advent of telescopes with large format imaging ar- 
rays, such as Herschel and JCMT, makes possible maps 
of star forming regions with large spatial dynamic ranges. 
It has become apparent in such images in both gas 
(Narayanan et al. 2008) and dust emission (Men'shchikov 
et al. 2010; Arzoumanian et al. 2011) that, when prop- 
erly sampled, molecular gas clouds exhibit filamentary 
structures on many spatial scales. This confirms the 
morphology that had been identified earlier in studies 
of visual extinction (Schneider & Elmegreen 1979). The 
origin of such elongated structures is unclear. A filament 
as a projected two dimensional structure, can be consis- 
tent with different theories. In particular, models suggest 
that gravitationally bound cigar-like clouds will collapse 
to become spindles, while pancake-like clouds will con- 
tinue to flatten under the influence of their self-gravity 
to become sheet-like structures (Lin et al. 1965; McKee 
& Ostriker 2007). In shock- induced molecular cloud for- 
mation models, the region around shock front will also 
be flattened (Vazquez-Semadeni et al. 2006). 

To distinguish a true filamentary (cylinder like) struc- 
ture from a sheet seen projected essentially edge-on, we 
have sought to obtain the line of sight dimension of a 
prominent dense filament in Taurus, B213 (Fig.[T]). Tau- 
rus is one of the closest star forming regions. The por- 
tion of Taurus molecular cloud that contains B213 fil- 
ament has been covered by a rich collection of surveys 
at optical through radio wavelengths (Goldsmith et al. 
2008; Rebufl et al. 2010). For this study, we have ob- 
tained data along "cuts" in the direction orthogonal to 
the elongation of the filament, in order to evaluate the 
variation of the excitation conditions in the filament, and 
to determine its volume density from which the line of 
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sight dimension can be determined. 

The column density of the molecular gas can be deter- 
mined by optically thin tracers including gas and dust. 
Goldsmith et al. (2008) presented a column density map 
of the whole Taurus region through analysis combining 
^■^CO 1-0 and ^^CO 1-0 transitions. There have been only 
a few previous attempts to determine the volume density 
of gas in Taurus. The first of these was by Avery et al. 
(1982), who observed a few positions towards TMC-1, 
a well-known high-density region. Avery et al. (1982) 
obtained five transitions (J=:4-3, 5-4, 9-8, 10-9, and 11- 
10) of cyanoacetylene (HC3N) at millimeter wavelengths. 
Their excitation analysis found a halo component having 
density 8x10'^ cm~^ and a core having a density 6x10 
cm~^. TMC-1 is part of the observed ring-like structure, 
but is not really a filament as are many other structures 
seen in Taurus. Schloerb, Snell & Young (1983) studied 
the same region in three transitions of IIC3N, and found 
a single density 10^ cm~^ to flt their data satisfacto- 
rily. If the '^l to 2 arcminute telescope beams are filled 
by gas at these relatively high densities, the line of sight 
dimension must be quite small, about O.lpc, in order to 
be consistent with the column density. Considering the 
molecular ring to be a tube of circular form, the diameter 
of the tube is no more than 0.4 pc. This does not give 
clear evidence for any extreme geometry. Pratap et al. 
(1997) derived the density for positions in TMCl based 
on the J=4-3, 10-9, and 12-11 transitions of HC3N. The 
range of densities are generally consistent with the core 
value from Avery et al. (1992), but Pratap et al. did not 
observe TMCIC. Onishi et al. (2002) identified a large 
number of cores in Taurus using H^^CO^ J = 1-0. The 
detection of this spectral line is generally suggestive of 
high densities, 10^ cm^'^, but since only a single transi- 
tion of this species was observed, the density could not 
be directly determined. Stepnik et al. (2003) made a 
cut across a dense filament in Taurus at multiple wave- 
lengths between 200 to 400 microns. Their model (which 
has large uncertainties) suggests a peak density ~ 6x10"^ 
cm~^, which suggests a true filamentary geometry. 
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Fig. 1. — ^^CO column density map of the Taurus molecular cloud containing B213 filament. The coordinates in this and subsequent 
figures are J(2000). The red letter "A" indicates the location of HC3N measurements and is the same as the position A in Fig. |4] The 
red square indicates TMCIC. 



In this study, we focus on using J=2-l and J=10-9 
transitions of HC3N as a density probe, which has the 
following advantages. The excitation of these spectral 
lines is sensitive to the volume density in likely range of 
this quantity of relevance to Taurus. When the 2-1 tran- 
sition at 18.2 GHz is observed with the 100 meter Green 
Bank Telescope (GBT), the resulting FWHM beam size 
is close to 40", which is similar to the 69" beam of the 12 
meter telescope of the Arizona Radio Observatory (ARO) 
at the 90.979 GHz frequency of the 10-9 transition. 

2. OBSERVATIONS 

At the Green Bank Telescope (GBT), we obtained data 
on the HC3N J=2-l transition toward positions across 
several filaments and TMCIC in the Taurus molecular 
cloud. The data were taken in two runs on 12 and 25 
February 2009. The observations were carried out in fre- 
quency switching mode with a frequency throw of 3 MHz, 
which is relatively small compared to the total bandwidth 
of 50 MHz. The total number of channels is 8096 with 
a velocity resolution of 0.1 kni/s. The resulting spectra 
have flat baselines and the emission is easily recovered 
through folding procedures (Fig. ^ . The intensity cali- 
bration was done following the procedure in the GBTidl 
manual written by J. Braatz. The calibration requires 
the knowledge of atmospheric zenith optical depth at 18 



GHz, which was provided by the GBT weather forecast 
tool to be 0.0183 for 12 February and 0.035 for 25 Febru- 
ary. The resulting antenna temperatures were further 
scaled to main beam temperatures by using a main beam 
efficiency of 0.88 (provided by GBTidl). 

At the Arizona Radio Observatory 12m telescope 
(ARO), we observed the HC3N 10-9 and 12-11 transi- 
tions toward the same set of filaments. Only the 10-9 
transition was detected and only in positions across the 
B213 filament. Between 2009 and 2010, we also carried 
out three runs to map out the spatial extent of the 10-9 
emission in B213. The HC3N data were taken in posi- 
tion switching mode at ARO. A simple correction was 
then applied to scale the data to main beam antenna 
temperature using the main beam efficiency of 0.87 at 91 
GHz from the ARO online manual. As seen in Fig. |31 
the hyperfine components of the 10-9 transitions are not 
resolved. 

3. RESULTS OF DATA ANALYSIS 

In all observing runs at both GBT and ARO, we ob- 
served TMCIC regularly to check the system and use 
it as a secondary flux calibrator. The analysis of these 
spectra are discussed in the following subsection using 
TMCIC as an example. The scientific analysis in SQland 
ij5]will focus on observations of B213. 
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3.1. TMCIC 

TMCIC exhibits strong HC3N J=2-l emission with 
peak main beam antenna temperature of the F=3-2 com- 
ponent being over 3 K. We did not detect the 2-1 hne 
toward other targets in Taurus except for the B213 fila- 
ment. This is a reflection of a much lower abundance of 
HC3N compared to that of TMCIC. The chemical differ- 
entiation in the Taurus region, although an interesting 
topic by itself, is beyond the scope of this paper. The 
dense core TMCIC has been known to exhibit strong 
emission in transitions of carbon chain molecules. We 
observed this source both as a relative intensity calibra- 
tor and as a test of our analysis procedures. 

The HC3N 2-1 transition has 6 hyperfine compo- 
nents, five of which can be clearly seen in our spec- 
trum of TMCIC. We fit all hyperfine components using 
the"hfs" fitting method in the GILD AS CLASS package 
(http://www.iram.fr/IRAMFR/GILDAS). The basic as- 
sumptions are a single line width for all hyperfine compo- 
nents and no blending, both of which are reasonable for 
TMCIC. The fitted line characteristics are given in Table 
2. The fitted peak optical depth is 0.47. This is a modest 
opacity, estimated with the reasonable assumption that 
excitation temperature of all hyperfine components are 
the same. 

The optical depth of the F=3-2 hyperfine component, 
and the corrected peak antenna temperature of the F=3- 
2 component are given in columns 7 and 5 of Table 2. 
The integrated intensity of each hyperfine component is 
corrected for its own optical depth, and the corrected 
integrated intensities for the J=2-l transition are given 
in column 6 of Table 2. 

3.2. The B213 Filament 

The data toward B213 were taken in the same setup 
as used for TMCIC. We observed a number of positions, 
which are indicated in Fig. 01 which shows all observed 
lines of sight toward B213 in both the J=2-l and the 
J=10-9 transition. We have detected both 10-9 and 2-1 
transitions in a total of 5 adjacent positions. The coordi- 
nates for these positions are given in Table 1. The extent 
of emission is delineated on three sides by non-detections 
and is well constrained within the ridge. There are two 
velocity components of HC3N in this region, which corre- 
spond to components seen in ^"^CO. The main component 
is at around 5.7 km/s, while the second component peaks 
at around 6.8 km/s. In the direction orthogonal to the 
elongation of the filament, the size of HC3N emission is 
between 2 and 3 arcminutes. The analysis of the density 
structure will focus on these 5 positions. 

These spectra are complicated by the fact that there 
are two velocity components in this filament. The dis- 
placement of these two components is about 1 km/s, 
which makes them blended in tracers with larger velocity 
dispersion, such as CO and its isotopologues. Due to the 
narrowness of HC3N lines, the two velocity components 
are clearly resolved. In Fig. [SJ we show the J=10-9 line 
at position A. A two component Gaussian fit reveals the 
peak velocities to be 5.6 and 6.7 km/s. The peak tem- 
peratures are 0.20 and 0.16 K, and the FWHM 0.28 and 
0.20 km/s, respectively. 

The J=2-l line of HC3N in B213 contains both hyper- 
fine and multiple velocity components. In Fig. [6] we show 
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Fig. 2.— Spectrum of the HC3N J=2-l line toward TMCIC. 
A double-frequency-switching mode was employed to obtain these 
data. The negative peaks (blue labels) are a consequence of the 
folding process used to analyze a frequency switched spectrum. 
The negative F=2-2, F=l-0, F=2-l, and F=3-2 components are 
a result of a shift to higher frequency, while the F=l-1 negative 
component is due to a shift to lower frequency. The intensity scale 
is the main beam antenna temperature, as described in the text. 
The hyperfine fit to the spectrum gives peak antenna temperature 
(corrected for optical depth) of the F=3-2 component to be 5.3 K, 
Vis,, equal to 5.4 km/s, and FWHM equal to 0.21 km/s. 
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Fig. 3.— Spectrum of the HC3N J=10-9 Une toward TMCIC. A 
single Gaussian fit gives Vjs,- equal to 5.2 km/s, FWHM equal to 
0.29 km/s, and integrated area equal to 0.38 K km/s. The peak of 
the fitted Gaussian is 1.23 K. 

the J=2-l spectrum of position A, in which a total of four 
peaks are visible due to the two stronger hyperfine com- 
ponents for each of the two velocity components, on the 
same velocity range as that of Fig. [S] The two observed 
hyperfine components of the 2-1 line are separated by 
about 1.5 km/s. We have performed hyperfine fits for 
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Fig. 4. — A summary of our pointings and HC3N detections. The 
underlying image is of magnitude of 2MASS extinction (Pineda et 
al. 2010). The blue symbols represent non-detections, the green 
symbols represent detection of the 5.6 km/s component, and the 
red symbols represent detection the 6.7 km/s component (see text). 
The circles denote the J=2-l transition and the squares denote the 
J=10-9 transition. The combination of shape and color represents 
whether emission is a particular velocity component is detected or 
not in a particular transition. 
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Fig. 5. — Spectrum of the HC3N J=10-9 line toward position A 
in the B213 filament. The red line is the result of a two velocity 
component Gaussian fit with parameters given in the text and in 
Table 2. 



each of the velocity components to derive their optical 
depths. The lower velocity component is found to have 
a peak optical depth of 0.48, similar to that of TMCIC. 
The higher velocity component is found to be optically 
thin. The total column densities of TMCIC and B213 
traced by 2MASS extinction (Pineda et al. 2010) are sim- 
ilar. The difference in line strength by a factor of 10 or 
more is thus attributable to a much higher HC3N abun- 



TABLE 1 
Source List 



Source 


Right Ascension 


Declination 




(J2000) 


(J2000) 


TMCIC 


04:41:38.82 


26:00:42.3 


B213 A 


04:18:38.45 


27:23:00.6 


B213 B 


04:18:35.45 


27:22:20.6 


B213 C 


04:18:38.45 


27:21:40.6 


B213 D 


04:18:41.45 


27:22:20.6 


B213 E 


04:18:38.45 


27:22:20.6 



dance in TMCIC. This chemical variation is not a focus 
of this paper and does not affect the excitation analysis 
and the derived volume density in the following sections. 
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Fig. 6. — Spectrum of the HC3N J=2-l line toward position 
A in the B213 filament. The velocity scale is referenced to the 
strongest (F=3-2) hyperfine component, which shows two velocity 
components, at 5.7 km/s and 6.8 km/s. The peaks at 7.3 and 
8.5 km/s are due to the F = 2-1 hyperfine component. The blue 
and red lines represent the hyperfine fits including two velocity 
components. 



4. DENSITY DETERMINATION THROUGH 
EXCITATION ANALYSIS OF HC3N 
TRANSITIONS 

Filaments can be defined as elongated structures in 
two spatial dimensions observed in images of both dust 
and gas. Such structures are found to be prevalent in 
star forming regions. The physics of filament formation 
and its implication for star formation hinges on the three 
dimensional geometry of such structures. The essentially 
unknown line of sight (los) dimension of structures that 
have very unequal dimensions in the plane of the sky 
(pos) is critical for determining whether the structures 
are filaments if djos is comparable to the smaller pos 
dimension or sheets viewed edge on if dios is comparable 
to the larger pos dimension. To determine djos , one need 
both volume density and column density. In this section, 
we derive the volume density based on excitation analysis 
of HC3N transitions. 



5 



In order to compare our data with statistical equilib- 
rium calculations of the excitation of the observed HC3N 
lines, we correct the integrated intensity of the J =2-1 line 
by the optical depth derived through fitting the hyperfine 
components. The integrated intensity of each hyperfine 
component is corrected for its own optical depth. The 
line ratio in Table 2 refers to the ratio between the total 
J=2-l integrated intensity and the integrated intensity of 
the J=10-9 line. This line ratio is what is modeled by the 
intensities given by our radiative transfer calculations in 
which the hyperfine components are not explicitly con- 
sidered. 

The situation for analyzing collisional excitation of 
cyanoacetylene in order to determine the local density 
has improved significantly of late. Green and Chapman 
(1978) carried out close-coupling calculations of rota- 
tional excitation in collisions with helium representing 
molecular hydrogen in its lowest, J = 0, spherically sym- 
metric state. These calculations were used in some of 
the works cited above. The first calculations that treated 
ortho- and para- H2 as distinct species were carried out 
by Wernli et al. (2007a). The improvements in the cal- 
culation of the potential energy surface were expected to 
yield improved accuracy for the collision rates. However, 
as explained by Wernli et al. (2007b), problems in the 
numerical realization resulted in the results for ortho-H2 
being invalid, while those for para-H2 are thought to be 
satisfactory. 

Quantum calculations by Faure and Wiesenfeld (2011) 
of collisions between ortho-H2 and HC3N have recently 
been made available to us. Collisions with para-H2 show 
a propensity for A J = 2 collisions, but the rates are only 
modestly smaller for A J = 1 and A J = 3. The new rates 
for collisions with ortho-H2 are quite different, favoring 
A J = 1. The result is that while the largest ortho-H2 
rates are up to a factor of 4 larger than the para-H2 rates, 
the overall collision rates are more nearly equal. Both the 
ortho-H2 and para-H2 rates are modestly larger than 
the Green and Chapman rates, even when the latter are 
scaled by a factor 1.39 to account for the lower reduced 
mass of H2 compared to that of He. 

We have used the RADEX code (Van der Tak et al. 
2007) to calculate the statistical equilibrium level pop- 
ulations and resulting line intensities from a molecular 
cloud. We adopt a column density N(HC3N) = IQi^ 
cm~^ and a line width of 1 km/s, which ensures that 
all transitions are optically thin, and which is close to 
A^(HC3N) determined for the B213 filament, in what fol- 
lows. In Fig. [7]we show the ratio of the intensity of the J 
= 10-9 to that of the J — 2-1 transition as a function of 
molecular hydrogen density, for a kinetic temperature of 
10 K. For optically thin emission and equal line widths 
(as found in the present data), the ratios of peak and 
integrated temperatures are the same. As seen in Fig. [T) 
the line ratio reaches specified value for a lower density 
of ortho-H2 than of para-H2. The difference, however, 
is only about a factor of 1.3 to 1.5 in density. 

The ortho to para H2 ratio (OPR) is quite uncertain. If 
the molecular hydrogen is formed with its internal levels 
populated at a temperature comparable to the molecule's 
binding energy, the OPR should just reflect the relative 
statistical weights of the two spin modifications, since 
their energy difference is so small in comparison. This 



would result in an OPR equal to 3. On the other hand, 
equilibration at the cloud temperature should result in a 
ratio determined by the J=l to J= energy difference. 
Since this is equivalent to 171 K, the result is an OPR 
close to zero for a dark cloud kinetic temperature of 10 K 
- 15 K. Observational results are difficult to obtain, and 
are mixed, but suggest that an OPR ~ 1 is a reasonable 
upper limit even in moderately warm regions (Neufeld, 
Melnick & Harwit 1998). This may be a reflection of 
timescales for interconversion of the spin modiflcations. 

Due to the very modest difference between the ortho 
and para collision rates, even though one might have to 
consider OPR between and 1, the excitation of HC3N is 
only very weakly dependent on this poorly-known quan- 
tity. In dark clouds there is evidence for enhanced molec- 
ular deutcration, which requires very low ortho H2 abun- 
dance (Flower, Pineau Des Forets & Walmsley 2006). 
Given the above considerations, we will adopt the para- 
II2 rates (i.e. an OPR = 0) for obtaining densities in the 
B213 filament. 

Our imprecise knowledge of the kinetic temperature 
results in an additional uncertainty. ^^CO data suggests 
a kinetic temperature between 12 K and 15 K, but this 
may apply more to the outer portion of the filament, 
since the lower transitions of this species are very opti- 
cally thick and do not really sample the interior of the 
fllament. It is reasonable that the the inner portions in a 
region such as this with a H2 column density correspond- 
ing to ~ 10 mag of visual extinction will be significantly 
cooler (Evans et al. 2001; Shirley et al. 2005). Thus, for 
the optically thin IIC3N emission, a kinetic temperature 
of 10 K seems appropriate. In Fig. |S]we show the effect 
of changing the kinetic temperature from 10 K to 15 K. 
The hydrogen density inferred for a given observed ratio 
is lower for the higher temperature, but in the range of 
values observed (< 0.5), the density difference does not 
exceed a factor of 2. 

Other than the uncertainty in the knowledge of colli- 
sion rates, the uncertainty in the derived density is dom- 
inated by the uncertainty in the derived optical depth 
for the J=2-l transition. The uncertainty in the fit- 
ted optical depth is affected by both the RMS noise, 
which is 0.09 K and 0.03 K per channel for TMCIC and 
B213, respectively, and the line strength. For TMCIC, 
the sigma for the derived optical depth of the main hy- 
perflne component is 0.23 compared to its best flt value 
0.47. This corresponds to a one sigma range of J=2-l cor- 
rected integrated intensities between 2.3 K km/s and 3.0 
K km/s. The RMS for the fltted integrated intensities 
of the J=10-9 transition is considerably smaller, about 
3%. The range of line ratios for TMCIC is thus between 
0.17 and 0.13. The resulting uncertainty in density would 
also depend on the shape of the excitation curve (Fig.[71). 
The one sigma range of derived densities for TMCIC are 
1.25 xlO'* cm~'^ to 1.04x10^ cm~^, which corresponds to 
about a 20% relative uncertainty for the derived density 
for TMCIC. For the B213 positions, the statistical un- 
certainties for derived optical depth are generally on the 
same order as the value of optical depth itself. For po- 
sition A, the range of derived densities is 1.9x10"* cm~'^ 
to 2.9x10"^ cm~^, which corresponds to about a ±20% 
uncertainty. This is representative of the statistical un- 
certainty in the density determined for B213. 

In summary, the largest source of uncertainty in the de- 
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Fig. 7. — Dependence of the ratio of HC3N J = 10-9 to J = 
2-1 intensities as a function of H2 density. The molecular column 
density is chosen so that all transitions are optically thin; in this 
case the ratios for the peak and the integrated temperatures are 
the same. The three sets of points correspond to the Green & 
Chapman (1978) rates, the Wernli et al. (2007a) rates for para- 
H2, and the Faure & Wiesenfeld (2011) rates for ortho-H2. The 
observed antenna temperature ratio of 0.4 corresponds to a H2 
density ~ 10'* cm"''. 
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Fig. 8.— Ratio of HC3N J = 10-9 to J = 2-1 transition as a 
function of the H2 density, for kinetic temperatures of 10 K and 15 
K. The collision rates used are those of Wernli et al. (2007a) for 
para-H2 (J = 0). The reduction in the derived density that results 
when the observed intensity ratio is less than 0.5 is relatively small, 
less than a factor of 2, in changing the kinetic temperature from 
10 K to 15 K. 

rived density arises from the imprecisely known kinetic 
temperature, followed by the uncertainty in the H2 or- 
tho to para ratio. The uncertainty caused by noise in 
the line measurement is in the range 20 to 40 percent. 
Given the reasonable assumptions that have been made 
for each of the items as discussed above, we feel that 
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Fig. 9. — The measured volume densities for the 5.6 km/s (in 
green) and 6.7 km/s (in red) components at positions A through 
E in the B213 Taurus filament (Table 1), for a kinetic temperature 
of 10 K and OPR = (see text for discussion). The units of the 



density are 10 



The blue circles refiect the GBT beam size. 



The peak of H^'^CO"*" emission is indicated by the white circle, with 
diameter reflecting the beam size of the Nobeyama telescope and 
center corresponding to core 10b listed in Table 1 of Onishi et al. 
(2002). The dark vectors denote the orientation of the magnetic 
field Ijased on optical absorption measurements (Chapman et al. 
2011). The underlying image is as in Fig. |4] 



we have a determination of the volume density through 
the measurements of HC3N line ratio that is accurate to 
within a factor of ±50% (factor of three). 

The density we derive (next to last column of Table 2) 
for TMCIC is about 1.1 x 10"* cm'^. For B213, the de- 
rived densities are a factor of 1.4 to 2 higher due to the 
larger observed line ratio, except for position E, which 
has slightly lower density than TMCIC. The derived den- 
sities for the five spatial positions are indicated in Fig. |9l 

We also derive the column density of HC3N using the 
observed intensity and the excitation conditions derived 
above. For an optically thin transition, the column den- 
sity of the upper level is simply proportional to its inte- 
grated intensity 



Tmbdv 



(1) 



where Tmb is the main beam antenna temperature as 
described in the data section. We use the integrated 
intensities of the optically thin J=10-9 transition to ob- 
tain N{J — 10). To obtain the total column density of 
HC3N, we utilize the level populations provided by the 
same Radex calculations that produced Fig. [7] for colli- 
sions with ortho-H2. The resulting column densities are 
listed in the last column of Table 2. The HC3N column 
density of B213 is within a factor of two of 10^^ cm~^, 
while that of TMCIC is one order of magnitude higher. 
This is consistent with the well-known fact that carbon 
chain molecules have significantly enhanced abundance 
in TMCIC. 
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5. DISCUSSION: THE 3D STRUCTURE OF B213 

The volume density measurements enable us to exam- 
ine the 3D structure of the B213 filament. We first con- 
sider the 2D dimension of the B213 filament projected 
on the sky. At the distance of Taurus, one minute of arc 
corresponds to about 0.04 pc. The shorter dimension of 
the B213 filament in the neighborhood of position A is 
about 4 arcminutes or 0.16 pc as seen in 2MASS extinc- 
tion. The filament is elongated by a factor of at least 10 
in the orthogonal direction. 

The median 2MASS extinction for these positions is 
about 6.3 mag, which corresponds to A''(H2) = 6 xlO^^ 
cm~^. This is close to that found from ^^CO (Pineda et 
al. 2010) assuming X{CO) = 2x10-"*, indicating that 
there is relatively little depletion of CO, which if present 
would imply an excessively large total column density. 
The modest depletion is consistent with the density we 
derive and a depletion time scale of a few x 10^ yr (Bergin 
& Tafalla 2007; Pineda et al. 2010), if the filament age 
is not much more than that. The average density found 
from our HC3N observations is about (1.8 ± 0.7) x 10^ 
cm~^. From the column density and the volume density, 
the line of sight dimension is found to be 0.12 pc, which is 
comparable to the shorter dimension of the filament pro- 
jected on the sky. The derived line of sight dimension is 
one order of magnitude shorter than the filament's elon- 
gated dimension in the plane of the sky, which effectively 
rules out B213 being an edge-on sheet. 

As shown in Fig. 21 we detected the higher rota- 
tional transition (J=10-9) of HC3N in 14 spatial posi- 
tions, which are distributed contiguously along the fila- 
ment. Detection of the lower transition was limited in 
only 5 spatial positions not because of insufficient exci- 
tation, but due to the difficulty of detecting the multiply 
hyperfine-split J = 2-1 transition. Fig. |9] shows density 
measurements at 5 spatial positions. The two velocity 
components present in these positions yield very similar 
densities. All measured densities are within a factor of 
2.5 of each other over a projected separation of approxi- 
mately 0.1 pc. The line widths at all observed positions 
are small and their peak velocities are similar. The vol- 
ume density measured for these two velocity components 
cannot result from an isolated "dense core" only seen 
high density tracers such as H^'^CO^. These densities 
are also unlikely to be the result of only a few very dense 
cores much smaller than the beam and resulting beam 
dilution of their emission. We seem to be probing lines 
of sight representative of the B213 filament as a whole. 
We cannot, however, rule out completely a highly clumpy 
cloud structure consisting of a large number of very small 
dense clumps and an extended low density component. 
Evaluating such substantially more complex structures 
will require observations at much higher spatial resolu- 
tion. 

If B213 is indeed a self-gravitating cylinder, an ex- 
tended low density 'halo' resulting simply from the re- 
quirement of hydrostatic equilibrium should be present. 
Chandrasekharr & Fermi (1953), Larson (1985), Inutsuka 
& Miyama (1992) and others have studied the configura- 
tion and evolution of infinite isothermal self-gravitating 
cylinders. The radial density profile of such cylinder 
turns out to be of the Plummer type for both the sta- 
ble and collapsing case with a 'flat' (approximately uni- 



form density) central part and a halo characterized by 
an effectively rapidly-dropping power-law density n ^ 
(l/r)"* (see Eq. 1 in Inutsuka & Miyama 1992). In such 
cases, the density measured by IIC3N will probe only the 
central 'flat' part, as the 'halo' cannot excite the higher 
rotational transitions. The dimension thus derived can- 
not rule out the existence of the extended low density 
halo, but will characterize the central core of the flla- 
ment, which could be considered to the represent the 
radial size of the filament. 

Cylindrical, i.e. filamentary, structures are of increas- 
ing interest at the present time due to their prevalence in 
recent surveys of the ISM, e.g. Goldsmith et al. (2008), 
Andre et al. (2010), and Molinari et al. (2010). The fil- 
aments found in these surveys differ by orders of magni- 
tude in terms of physical length, mass, and their star for- 
mation content. These results show that, at least in 2D 
projection, elongated structures are a common feature of 
the interstellar medium. Filaments can be produced in 
numerical simulations of vastly different physical condi- 
tions. Filaments show up in a strongly magnetized (sub- 
Alfvenic) MHD simulation with ambipolar diffusion (Li 
et al. 2012). Filaments are also seen in simulation of su- 
personic converging flows without magnetic field (Gong 
& Ostriker 2011). Filamentary structures are not re- 
stricted to fluid simulations (with or without magnetic 
field, with or without dominant turbulence); filamentary 
structures are prevalent in cold dark matter cosmolog- 
ical simulations (Springel et al. 2005), where the most 
important physics are gravity and structure growth. 

To determine the driving factor of structure formation 
in the star-forming ISM, the measurement of physical 
parameters in addition to morphology is critical. The 
data presented here are a limited first step. One partic- 
ular relevant regime of ISM evolution is molecular cloud 
formation induced by colliding HI streams (Ballesteros- 
Paredes et al. 1999). A flattened structure (sheet) is 
expected at the shock front (Vazquez-Semadeni et al. 
2006). The subsequent structure growth depends on the 
relative importance of different unstable modes, which 
in turn depend on the physical conditions. In an illus- 
trative discussion by Heitsch et al. (2005), the structures 
that become visible largely reflect the ratio between vol- 
ume density and converging velocity. When the den- 
sity is high, thermal instability dominates and the dense 
structures that grow tend to be contained near the shock 
front, maintaining a fllamentary appearance. When the 
ratio of density to colliding velocity decreases, Kelvin- 
Helmholtz and nonlinear thin sheet instabilities grow 
faster and tend to produce structures perpendicular to 
the elongation of the shock front. 

We are not yet at a stage to differentiate among dif- 
ferent simulations, for which the methods themselves 
are evolving rapidly. It is, however, tantalizing to de- 
scribe the evolution of ISM leading to star formation 
as the following sequence: sheet — >■ filament — >■ cores. 
A flattened structure (sheet) can be produced by col- 
liding streams or other shocks as discussed above. A 
self-gravitating sheet will preferentially fragment into fll- 
aments (Miyama, Narita & Hayashi 1987). Inutsuka 
and Miyama (1992) carried out self-similar calculation of 
isothermal cylinder and derived a critical fllament "line 
mass", consistent with that in Ostriker (1964), beyond 
which the fllament is expected to fragment and collapse 
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into cores. This last stage of filament -to-cores evolu- 
tion seems to be emerging as a phenomenon often found 
in recent observational studies, such as Hacar & Tafalla 
(2011) and Schmalzl et al. (2010). 

In observing multiple transitions of HC3N toward Tau- 
rus filaments, we do not see any suggestions of a sheet 
extended along the line of sight. If a sheet had been 
present, it must have already evolved or dissipated. The 
line width and velocity gradient of IIC3N gas are both 
small at around 0.3 km/s. At the ~ 10 K temperature of 
gas in the filament, this is marginally supersonic. For a 
sound speed of 0.127 km/s corresponding to a line width 
of 0.3 km/s FWHM, the critical line mass is about 7.6 
Mq/pc. For B213, this corresponds to a column den- 
sity of N(H2) = 4x 10^^ cm~^, somewhat smaller than 
the observed value, suggesting that B213 is susceptible 
to collapse/fragmentation to produce cores. 

The line width of HC3N in B213 filament is similar to 
that of I80 in L1517 (Hacar & Tafalla 2011), while larger 
than that of N2H+ that is observed in the same study. 
Hacar & Tafalla also found a subsonic velocity difference 
among N2H"'" cores contained in the L1517 filaments and 
were able to fit the radial profiles of I80 emission using 
isothermal cylinders. These authors postulate a two step 
core formation processes including, first, a subsonic fil- 
ament, and second, fragmentation-produced cores. The 
measured volume density and line of sight dimension of 
B213 are consistent with it being a cylinder-like 'true' fil- 
ament in which HC3N traces the gas in a manner similar 
to I80 in L1517. The internal motions of the B213 fila- 
ment are still slightly siipersonic, which could mean that 
B213 is at an earlier evolutionary stage than L1517, hav- 
ing yet to go through the stage in which the line width 
drops and core formation begins. 

The volume density measurements and the velocity in- 
formation presented in this paper will be of additional 
value when radiative transfer is added to various simula- 
tions and modeled spectral line profiles can be compared 



directly with those observed. The density measurements 
in this study fit into a developing picture of molecular 
cloud evolution from turbulent gas to coherent filaments 
to subsonic cores. 

6. CONCLUSIONS 

We have determined the volume H2 densities along five 
lines of sight (which include a total 6 distinct compo- 
nents) in a prominent filament in Taurus, B213, to be on 
average (1.8 ± 0.7) x lO"' cm^''. The density determina- 
tions are based on excitation calculations that fit the line 
ratio between HC3N 10-9 and HC3N 2-1 transitions. The 
clearly-resolved hyperfine components of the 2-1 transi- 
tion enable us to determine the peak opacity of the 2-1 
line to be smaller than 1.3. The line of sight dimension 
of the dense portion of B213 filament is ~ 0.25 pc, com- 
parable to the smaller plane of the-sky dimension of the 
filament, which is about 0.16 pc. These dimensions are 
both much smaller than the longer dimension of B213, 
which is a few pc. The FWHM line wdith of HC3N is 
about 0.3 km/s, which is slightly supersonic. We con- 
clude that B213 is likely a true filament, but that it has 
not yet formed dense cores throughout. 
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TABLE 2 
Fitted HC3N Line Characteristics 



Source Comp. HC3N 2-1 HC3N 10-9 Line Ratio n(H2) N(HC3N) 







Vo 
km/s 


FWHM 
km/s 


T 

main 

K 


jTdv 
K km/s 


Tmairt 


Vo 
km/s 


FWHM 
km/s 


T 

K 


JTdv 
K km/s 


/ T{10-9)dv 
JT{2-I)dv 


10" cm3 


IOI2 cm-- 


TMCIC 


s 


5.4 


0.21 


5.34 


2.6 


0.47 


5.2 


0.29 


1.23 


0.37 


0.15 


1.1 


18. 


B213 A 


b 


5.7 


0.14 


0.45 


0.15 


0.48 


5.6 


0.28 


0.20 


0.058 


0.39 


2.3 


1.5 


B213 A 


r 


6.8 


0.17 


0.22 


0.083 


0.00 


6.7 


0.18 


0.16 


0.030 


0.36 


2.2 


0.80 


B213 B 


b 


5.7 


0.14 


0.14 


0.046 


0.00 


NA 


NA 


NA 


NA 


NA 


NA 


NA 


B213 B 


r 


6.8 


0.17 


0.20 


0.076 


0.00 


6.7 


0.19 


0.14 


0.029 


0.38 


2.3 


0.74 


B213 C 


s 


5.6 


0.17 


0.37 


0.14 


0.27 


5.5 


0.29 


0.10 


0.031 


0.22 


1.5 


1.1 


B213 D 


s 


5.6 


0.19 


0.24 


0.10 


0.14 


5.5 


0.19 


0.12 


0.024 


0.23 


1.6 


0.85 


B213 E 


s 


5.6 


0.17 


1.03 


0.39 


1.30 


5.5 


0.18 


0.19 


0.038 


0.097 


0.94 


1.9 



Note. — The units for the fitted FWHM hnewidth and central velocity Vo are km/s. The units for Tmain, the main beam antenna temperature 
from the fitting procedures described in the text, are Kelvin. In the second column, s denotes a single velocity component, b and r denote blue and 
red velocity components, respectively. Column 7 gives the peak optical depth of the F=3-2 (main) HPS component of the J=2-l transition. The 
integrated main beam temperatures having units of K km/s have been corrected for the optical depth where appropriate. An optical depth of zero 
indicates an optically thin line. "NA" indicates a nondetection. 



